Two-electron dephasing in single Si and GaAs quantum dots 
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We study the dephasing of two-electron states in a single quantum dot in both GaAs and Si. 
We investigate dephasing induced by electron-phonon coupling and by charge noise analytically for 
pure orbital excitations in GaAs and Si, as well as for pure valley excitations in Si. In GaAs, polar 
optical phonons give rise to the most important contribution, leading to a typical dephasing rate 
of ~ 5.9 GHz. For Si, intervalley optical phonons lead to a typical dephasing rate of ~ 140 kHz 
for orbital excitations and ~ 1.1 MHz for valley excitations. For harmonic, disorder- free quantum 
dots, charge noise is highly suppressed for both orbital and valley excitations, since neither has an 
appreciable dipole moment to couple to electric field variations from charge fiuctuators. However, 
both anharmonicity and disorder break the symmetry of the system, which can lead to increased 
dipole moments and therefore faster dephasing rates. 
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I. INTRODUCTION 

Both Si and GaAs quantum dot technologies are now 
weh-estabhshed as candidates for the implementation of 
scalable quantum computation.^ ^ Working qubits have 
been fabricated with several different architectures in- 
volving singly-occupied dots, including the single-spin 
qubit architecture pM^ the singl et-triplet scheme,^ ^ ^'^ 
and three-dot logical qubits!^^'^ 

Recently, we proposed a new "hybrid" qubit 
architecture^ with three electrons in two quantum dots 
that is potentially advantageous because fast qubit op- 
erations can be performed in a relatively simple archi- 
tecture. The hybrid design is capable of all-electrical, 
fast qubit gates, but the qubit wave functions have some 
charge character that gives rise to decoherence mecha- 
nisms that are not present in pure spin qubits. This 
means that characterizing decoherence becomes a most 
pressing issue, because charge decoherence is typically 
much faster than spin decoherence. 

In the singlet-triplet qubit, when there is a finite 
exchange-induced energy splitting between the singlet 
and triplet states, the dominant sources of decoherence 
have been found to be charge noise^^ ^^ and electron- 
phonon coupling. ^^ Much of the physics of the hybrid 
qubit is similar to that of the singlet-triplet qubit, with 
the main differences in the decoherence properties aris- 
ing because one of the dots contains two electrons. In 
the doubly-occupied dot, both electrons in a singlet can 
occupy the orbital ground state, while at least one of the 
triplet electrons must lie in an excited state. A further 
complicating factor in silicon-based quantum dot devices 
is the presence of two nearly degenerate, low- lying valley 
states .I^^H^ These levels are split near a sharp interface 
by an amount that is typically comparable to the orbital 
energy spacing. Hence, single-electron first excited states 
have two characteristic types: orbital, where the electron 
occupies the same valley state but the P-like first ex- 



cited state of the lateral confinement potential, and val- 
ley, where the electron is in the orbital ground state and 
a higher valley state. 

In this paper, we calculate singlet-triplet dephasing 
rates in a doubly-occupied quantum dot, extending the 
spin relaxation calculations previously done for both 
GaA&^ and Si.-^ Unlike relaxation, pure dephasing is 
due to processes that conserve spin and does not in- 
volve energy exchange with the environment. For both 
electron-phonon coupling and charge noise, we consider 
the limiting cases of purely orbital (for both GaAs and 
Si) and purely valley (for Si) excited states. 

We find that for GaAs, polar optical phonons are the 
main source of dephasing, leading to a decoherence rate 
of ~ 5.9 GHz. For Si, the phonon-mediated dephasing 
rate depends on the type of excitations supported by the 
quantum dot. For a first excited state that is orbital-like, 
intervalley optical phonons lead to a decoherence rate of 
~ 140 kHz. For a valley-like first excited state, this same 
phonon channel results in a faster decoherence rate of 
~ 1.1 MHz. For a perfectly harmonic, disorder- free dot, 
we find that dephasing due to charge noise is strongly 
suppressed. This is because the effective dipole moment 
between the singlet and triplet states vanishes for both 
orbital- and valley-like excitations. If we allow for anhar- 
monicity and an effective dipole moment, we find that 
both phonon and charge noise dephasing channels are 
of similar strengths in Si, but phonons are the limiting 
mechanism in GaAs. Assuming a gate operation speed of 
10 GHz (quite feasible for the hybrid qubit), the decoher- 
ence rate in silicon is consistent with the achievement of 
10^ operations per coherence time, while the decoherence 
rate in GaAs is too fast for viable hybrid qubit operation. 

This paper is organized as follows. In Sec.jllj we briefiy 
review the quantum states that are relevant to the system 



considered in this paper. Next, in Sec. |III|we formulate 
the problem of intra-dot singlet-triplet dephasing due to 
the electron-phonon coupling, following the formalism of 



Ref. [191 We consider first GaAs, then both pure orbital 
and valley excitations in Si. In Sec. \LV\ we calculate 
dephasing due to charge noise, and compare to phonon- 
induced dephasing. Finally, in Sec.[V|we discuss the role 
that the dephasing mechanisms we have addressed are 
likely to play for qubits and suggest methods for miti- 
gating decoherence. 



II. TWO-ELECTRON STATES IN A QUANTUM 

DOT 

The lowest energy eigenstates of two electrons in a sin- 
gle quantum dot (neglecting spin-orbit interaction) are 
the singlet 



\S) = iVo) 
and the three triplet states 
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where IV^o) is the two-electron spatial ground state and 
IV^i) is the first excited state. We work within the Heitler- 
London approximationf^ in this approximation, |V^o) ^ 
1 00), indicating that both electrons are in their single- 
electron ground states, and |?/^i) ^ (|01) — |10)) />/2, 
where one electron is in its ground state and the other is 
in its first excited state. This paper considers dephasing 
due to the elect ron-phonon coupling and charge noise, 
both of which conserve spin, so from here on we will fo- 
cus on the spatial component, with the appropriate spin 
wavefunction understood. We refer to the triplets col- 
lectively as \T) when the particular spin configuration is 
not important. 

In GaAs, the wavefunctions of the ground and first ex- 
cited states have identical dependence on z (the direction 
perpendicular to the quantum well), but have S- and P- 
like transverse envelopes in the x-y plane (the plane of 
the quantum well). For a quadratic quantum dot con- 
finement potential, the effective mass approximation for 
the ground state wavefunction is 

^o^°^^(r) = «(r)^e-(-'+^')/(2^') (3) 

VttL 

while the first excited state wavefunction is 

^GaA.(r) = „(r)ffl(f±Me-(^^+.^)/(2i^^), (4) 

where u{r) is the periodic (with the lattice periodicity) 
component of the Bloch function at the conduction band 
minimum (the F-point in GaAs), F{z) is the envelope 
function along z, and L is the lateral extent of the wave- 
function. 



For silicon, we write the ground state wavefunction as 

^o^^(r) = ^_(r)^e-ii-(-^+^^), (5) 



where 
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and u± (r) is the periodic part of the Bloch function eval- 
uated at the conduction band minimum located at ±/coi. 
In a quantum dot fabricated in a strained silicon quan- 
tum well, ko ~ 0.82 x 27r/a, with a = 0.543 nm the length 
of the Si cubic unit cell.^^ Depending on the magnitude 
of the valley splitting introduced by the sharp interfaces 
and the electric field in the z-direction, the lowest energy 
excited states can either be valley-like or orbital-like. For 
the case of an orbital-like excitation (large valley split- 
ting), the wavefunction is 
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When the valley splitting is smaller than the or- 
bital splitting, both the ground state and first excited 
state have S-like transverse wavefunctions, but their z- 
direction wavefunctions are different valley states. The 
ground state wavefunction is still given by Eq. ([5|, but 
the first excited state is now 



^1 (r) 
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The periodic parts of the Bloch functions, 'u(r), have 
discrete Fourier spectra, with contributions occurring at 
reciprocal lattice vectors G. Hence, when performing 
calculations with the full wavefunctions, as we do in 
this paper, it is convenient to decompose u{y) into the 

g^j^J26l27l3Qn 
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where the expansion coefficients a(G) are independent of 
r. For low- frequency dephasing channels such as acous- 
tic phonons and charge noise only the G = mode con- 
tributes significantly.'^ For high-frequency processes such 
as optical phonon couplings, contributions with G ^ 
can also be important. We assume that the optical 
electron-phonon couplings are independent of G. With 
this assumption, one can prove that the calculation of de- 
phasing rates is independent of the form of the periodic 
part of the Bloch functions. Hence, for all instances we 
consider in this paper, we may ignore the periodic part 
of the Bloch functions. 

In a real system, disorder would cause the excited 
states to have mixed valley and orbital characteristics,^ 
for which pure orbital- like and pure valley- like first ex- 
cited states represent limiting cases. Hence, it is impor- 
tant to consider both the pure valley and orbital excita- 
tions described above. 



III. DEPHASING VIA THE 
ELECTRON-PHONON INTERACTION 

We now consider the dephasing of two-electron states 
in a single quantum dot due to the electron-phonon inter- 
action, following the techniques of Ref. [19l In Sec. |III A[ 
we consider dephasing in GaAs due to deformation poten- 
tial, longitudinal and transverse piezoelectric, and polar 
optical phonons. In Sec. |IIIB[ we turn to silicon, where 
we may have either valley or orbital excitations, and the 
relevant dephasing channels are through intravalley de- 
formation potential and intervalley optical phonons. 

The general form of the electron-phonon interaction 
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where each G is a reciprocal lattice vector, q is con- 
strained to the first Brillouin zone, p is the electron den- 
sity operator, a^^x and a^ ^ annihilate and create, respec- 
tively, phonons with wave vector q, and A indexes the 
phonon mode. Since we will treat the the relevant modes 
separately, we will suppress this sum over phonon modes 
in the calculations that follow. The electron-phonon cou- 
pling, Mq+G, is defined by 



M(q+G) = -K»(q+G) [(q+ G) • C] 



2pmfi^q 



(11) 



where Vd is the electron- ion potential, ^ is the phonon 
polarization vector, pm is the crystal mass density per 
unit volume, Q is the crystal volume, and cjq is the 
phonon frequency. Since the singlet and triplet have dif- 
ferent charge distributions, they are dressed differently 
by the phonons. The phonons can themselves decohere, 
which in turn causes dephasing between the singlet and 
triplet states. Following Ref. 19, the singlet-triplet de- 
phasing rate due to the electron-phonon coupling is 
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(12) 
where 7q is the phonon relaxation rate, and A is a 
Fourier component of the charge density difference be- 
tween triplet and singlet states: 

^(q +G) = \ m Kq + G) IT) - {S\ Kq + G) \S)) . 

(13) 
Using the approximate forms for the singlet and triplet 
states detailed in Eqs. (fTl) and (pi), we have 

A(q + G) = 1 ((1| e'('i+G)-'- |1) - {0| e^('i+«)-'- |0)) , 

(14) 
where |1) is the single-particle first excited state and |0) 
is the single-particle ground state. In the following sub- 
sections, we evaluate Eq. (12) for the different types of 
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FIG. 1: (Color online) Calculated singlet-triplet dephasing 
rates of two electrons in a single dot via the phonon mech- 
anisms considered in Sec . |III| and the charge noise mech- 
anism considered in Sec. |IV[ (a): Plot of dephasing rates 
versus the lateral extent of the electron wavefunctions for 
the mechanisms pertinent for GaAs: polar optical phonons 
(Eq. (24)), deformation potenti al ph onons (Eq. (18)), trans- 
verse piezoelectric phonons (Eq. ([22|), and longitudinal piezo- 
electric phonons (Eq. (20)). (b): Plot of singlet triplet de- 



phonons in both GaAs and Si. 



phasing rates versus lateral wavefunction extent for different 
mechanisms in Si: intervalley phonons for valley excitations 
(Eq. (33)), and interv alley phonons (Eq. ( |32| )), longitudinal 
acoustic phonons (Eq. (29)), and transverse acoustic phonons 
(Eq. (30)) for orbital excitations. The dephasing rates due 
to charge noise plotted for both materials systems are deter- 
mined from Eq. ( |40[ ) , using the ene rgy fluctuations for orbital 
excitations in a single dot (Eq. (36 )), and assuming zero dipole 
moment. 



A. Phonon-induced dephasing in GaAs 



For the purely orbital excitations supported by GaAs, 
the ground state transverse wavefunction is an S-orbital, 
while the first excited state transverse wavefunction is a 
P-orbital, as given in Eqs. ([3| and (p|. The three types 
of phonon couplings that contribute to decoherence are: 
deformation potential, piezoelectric, and polar optical. 
These differ only in the form of the electron-phonon cou- 
pling M, so their calculations proceed similarly. 

In all cases, we assume a Gaussian form for the wave- 



function in the z-direction: 
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where d is the confinement length along the growth axis, 
which is typically a few nanometers. Choosing this form 
for the wave function represents an approximation, but 
it captures the relevant physics and allows us to obtain 
analytic results. Using this approximation and Eqs. O) 
and Q, we use Eq. ([l4| to obtain 
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(16) 
Here, we have taken G = 0, since 1^41 goes to zero 
rapidly for |q| = q ^ l/<i, 1/i^, both of which are much 
smaller than the size of the first Brillouin zone. 

We first consider deformation potential phonons, for 
which the electron-phonon coupling M has the fornJ^ 



M^,^^M=Dq 
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where D = 8.6 eV is the deformation potential con- 
stant and pm = 5.33 x 10^ kg/m^ is the mass density 
of GaAs. The angular frequency cjq is given by the stan- 
dard relationship for acoustic phonons: cjq = VgQ^ where 
Vs = 5.2 X 10^ m/s is the longitudinal speed of sound in 
GaAs, averaged over direction. ^^ 

The last piece of information we need to compute 
the dephasing rate is the phonon relaxation rate 7^. 
Eq. (16) implies that \A\'^ is strongly peaked, so it de- 



termines the q-values that contribute to the integral in 
Eq. (12). In GaAs, |Ap goes to zero both as q goes to 



zero and when q ^ l/d, 1/L. Hence, the low- and high- 
frequency behaviors of 7^ are not important. For suffi- 
ciently high frequencies at low temperature, it is expected 
that two frequency- dependent phonon attenuation chan- 
nels will become relevant: anharmonic decay and isotope 
scattering. ^^ However, the frequencies we consider here 
are low enough that these mechanisms are unimportant, 
and the dominant source of phonon relaxation is due to 
interface scattering.^ 

To obtain an estimate for the phonon relaxation rate, 
we use experimental measurements of phonon attenua- 
tion due to interface scattering, which were performed 
at low-temperatures in Si.'^ To convert between the 
two, one uses 7^ = 2aqVsr^ where aq is attenuation 
and Vs = 8.49 x 10^ cm/s is the speed of sound in Si 
along [100]. For LA phonons in Si, low-temperature 
measurements have shown that at low frequencies (up 
to 100 GHz), phonon attenuation is roughly frequency- 
independent, and is about 2.5 cm~^ along the [100] 
direction. ^^ This translates to a low-frequency experi- 
mental limit of Jq"^ = 4.25 MHz, which will serve to 
give us an estimate on the phonon relaxation rate. Since 
this mechanism is due to the geometry rather than the 
particular material properties of Si, we will also use the 
above relaxation time for acoustic phonons in GaAs. 



By switching to polar-cylindrical coordinates, evalua- 
tion of the integral in Eq. (12) is straightforward. We set 
(^^ + (7^/2)^ ^ cjq, which is valid because the frequencies 
that contribute to the integral satisfy cjq ^ 70. Electro- 
statically defined quantum dots typically obey L ^ d^ 
so we expand the integration result to first order in d/L^ 
obtaining an expression for F^J ^' , the singlet-triplet 
decoherence rate in GaAs due to deformation potential 
electron-phonon coupling: 

GaAs,DP ^ D\^H2L/d)-Z)^^^ _ 6.9 2 
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Here, the dependence on L and d can be understood by 
power counting in Eq. (12). For the polar couplings we 



consider next, the L and d dependencies can be under- 
stood by examining the qz integral in Eq. (12) over the 
range where \qz\ <C \qx\ , |%|, followed by power counting. 
We next consider piezoelectric coupling, which con- 
tribute in both longitudinal and transverse phonon 
modes. In this case, the electron-phonon coupling is^ 
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where e is elementary charge and 614 = 1.38 x 10^ V/m 
is an elasticity tensor component. For longitudinal 
phonons, ^ = q/g. Integration of Eq. (p^, ex- 



panded to lowest order in d/L, yields F^^ *' ' , the 
singlet-triplet dephasing rate in GaAs due to piezoelec- 
tric coupling between electrons and longitudinal acoustic 
phonons: 
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(20) 



In the limit when the media is considered isotropic and 
homogeneous, the two transverse phonon branches are 
degenerate, and we can choose any two orthogonal po- 
larizations. One possible polarization is 
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and any rotation of this vector about q is also a valid 
transverse polarization. We average over this plane be- 
fore integrating over q. This complicates the resulting in- 
tegral, but it can still be carried out analytically, yielding 
the singlet-triplet dephasing rate in GaAs 
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due to piezoelectric coupling between electrons and trans- 
verse acoustic phonons: 
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For polar optical phonons, the electron-phonon cou- 
pling is^ 



M^fA.(q) 
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where eoo = 10.89e^ac and eo = 12.9eyac are the high 
and low frequency limits of the GaAs dielectric func- 
tion, and Cyac is the vacuum permittivity. The frequen- 
cies of optical phonons are essentially g— independent, 
with hwo = 36.35 meV. Since optical phonons are much 
higher in frequency than acoustic phonons, they also 
have much shorter lifetimes, with measurements indicat- 
ing Jq^ ^ 160 GHz.^^ The integration proceeds similarly 
to the acoustic cases, and the resulting singlet-triplet 
dephasing rate T^^ ^' in GaAs due to polar optical 
electron-phonon coupling is, to first order in d/L: 
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In Fig. [T] we plot the four dephasing rates considered 
in this section. Typical values and scalings are listed in 
Table [ij Polar optical phonons are the largest contribu- 
tion to dephasing, exceeding the others by at least five 
orders of magnitude. This is mainly due to the extremely 
fast decay of the high-frequency, optical phonons. 



B. Phonon- induced dephasing in silicon 

Unlike GaAs, Si quantum dots can support both val- 
ley and orbital excited electron states. The ground and 
first excited states in the case of an orbital excitation are 
given in Eqs. ([5| and ^. As we did for GaAs, we take 
the envelope in z to be Gaussian with width d. Then, 
defining Q = q + G, we evaluate Eq. (14) for orbital 
excitations in Si (^o*), obtaining 



^ ^-d\2ko+Ql)/2 ^ ^-d^i2ko-Ql)/2\ ^ (^25) 

where we used the fact that the three Gaussians are well- 
separated in Qz to drop cross-terms. If we instead have 
excitations as given by Eqs. ^ and (|8|, we evaluate 
Eq. (|14| for valley excited states in Si (Ay)^ obtaining 
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The expressions for \A^\^ and \A^''\^ in Eqs. ([25) 
and ([26]) above select the types of phonons that con- 
tribute significantly to dephasing through the integral 
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FIG. 2: (Color online) Cartoon of the absolute value of the 
Fourier transform of the difference between the triplet and 
singlet charge distributions, \A\ (Eq. ( [l3| ) in Si. Here, the 
shaded regions indicate the k— values that contribute signif- 
icantly to dephasing. Phonons that couple electrons in the 
same valley (intravalley processes) lie near the origin, while 
phonons that couple electrons in different valleys (intervalley 
processes) lie near k = ±2kokz. (a): The contribution re- 



sulting from an orbital-like first excited state (Eq. (25)). (b): 
The contribution resulting from a valley-like first excited state 
(Eq. ([26])). 



in Eq. (12). For orbital excitations, both phonons that 



couple electrons within the same valley and across val- 
leys contribute, and |Aq*P contains contributions from 
three toroids, each with peak radius a/2/L, situated in 
the kx — ky planes centered at Q^ = and Qz = ±2A:o. 
The phonons at the Qz = toroid correspond to in- 
travalley processes, where G = 0. The remaining two 
toroids at Qz = ±2/co are intervalley processes, where 
G = ^(47r/a)i. The regions of k— space relevant to 
orbital excitations are shown in Fig. [2] (a). For valley 
excitations, |Ay p contains contributions from two ellip- 
soids centered at Q = (0,0, ±2/co). Since the valley-like 
first excited state has the same envelope function as the 
ground state, long wavelength phonons cannot contribute 
to the singlet-triplet dephasing. This is clearly illustrated 
by the vanishing of |Ay p at small k (or long wavelength). 
The regions of k— space relevant to valley excitations are 
shown in Fig. [5] (b). 

Now that we have identified the most important 
phonon wave vectors for the different dephasing mecha- 
nisms, we discuss which electron-phonon coupling mech- 
anisms are most relevant. Since Si is not polar, the 
deformation potential electron-acoustic phonon coupling 
is the main contribution near the zone center (G = 0, 
q <C 27r/a). This coupling connects electrons to both 
longitudinal acoustic phonons, with matrix elemenlP^ 



\Mla{ci) 
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TABLE I: Typical values for the dephasing rates of the phonon-induced channels discussed in Sec. |IIH assuming a lateral 
electron confinement of L = 40 nm and a vertical confinement along the growth direction oi d — 3 nm. The channels depend 
on material properties and symmetry polar phonon couplings are absent in Si, transverse phonons do not couple electrons via 
the deformation potential in GaAs™ at low temperatures and biases GaAs has only one band minimum that participates in 
conduction, and valley excitations in Si are not connected by low-frequency phonons. The scaling column describes the primary 
dependence of the dephasing rate on L, the lateral wavefunction extent and d, the vertical wavefunction extent. In the scalings, 
(*) indicates a neglected logarithmic correction, while (**) indicates that the result applies in the limit of d/L -^ 0. 



Coupling mechanism 



Typical dephasing rate (Hz) 



GaAs Si: orbital excitations Si: valley excitations 



Scaling 



LA phonons 

Deformation potential 4.3 x 10^ 
Piezoelectric 6.0 x 10^ 



8.8 X 10^ 



TA phonons 

Deformation potential 
Piezoelectric 2.2 x 10^ 



5.2 X 10^ 



LO phonons 
Polar optical 
Intervalley 



5.9 X 10^ 



1.4 X 10"^ 



1.1 X 10^ 
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and to transverse acoustic phonons, with coupling^^ 
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where E^ = 5.0 eV and Eu = 8.77 eV are silicon deforma- 
tion potentials. As in Sec. |III A] above, for these acoustic 
modes we take the phonon relaxation rate, believed to 
be due to interface scattering, to be 70 = 4.25 MHz.'^ 
The intravalley piece of the orbital excitation is found 
by integration of Eq. (12), which proceeds very similarly 



to the GaAs deformation potential case we considered in 
Sec. |III A[ For deformation potential coupling between 
electrons and longitudinal phonons in Si, the singlet- 
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reciprocal lattice vectors that contribute significantly to 
the relevant integrals are G = ±{4:7r/a)kz^ which give 
q^ ^ T0.36(27r/a) = t4.17x 10^ m-^ The phonons that 
are responsible for this transition in silicon are due to g- 
type Umklapp processes. ^^ Although symmetry restricts 
these to be longitudinal optical phonons, experiments 
indicate that both transverse and longitudinal acoustic 
phonons participate through processes in which M(q) is 
first-order in q.^^ ^^ However, the acoustic phonons do 
not play a significant role here, both because their de- 
formation potential coupling to the electrons is weaker^ 
and they are much longer lived^^ than optical phonons. 
The LO phonons in Si have a nearly constant energy 
hwci = 62 meV.^^ The electron-phonon coupling arises 



triplet dephasing rate F^^ , to first order in d/L, is: f^^^ ^^ ^^^-^^^^ deformation potentialPEH 
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As was the case for GaAs, the dependence on L and d 
can be understood by power counting in Eq. (12). 



We perform the same averaging procedure as was done 
for the transverse phonons in the previous section and 
obtain F^^ , the singlet-triplet dephasing rate in Si due 
to deformation potential coupling between electrons and 
transverse acoustic phonons: 



|M(Q) 



LO\ 






(31) 



where the intervalley deformation potential Dif = 11.0 x 
10^ eV/cm.l^ 

Finally, we estimate the relaxation rate 7^ for opti- 
cal phonons in Si. As for optical phonons in GaAs, the 
short- wavelength longitudinal optical phonons that cause 
intervalley coupling have a much shorter lifetime than the 
long wavelength acoustic phonons that are responsible for 
intravalley coupling. The literature value we use for the 
relaxation rate is Jq^ = 118 GHz.!^ For the intervalley 
component of the orbital excitation, we get 
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— r kHz nm . 



(30) 



We next consider the intervalley contributions, which 
occur at kz = ±2/co, outside the first Brillouin zone. The 
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(32) 



Similarly, for the case of valley excitations we have 
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(33) 



Fig. [T] shows the dephasing rates for both GaAs and 
Si. Typical values and scalings with L and d are listed 
in Table [Tj In Si, as in GaAs, most of the dephasing is 
due to the fast decay of optical phonons: in silicon these 
high-frequency phonons couple electrons across valleys. 

So far, we have only considered pure valley and pure 
orbital excitations. However, for non-ideal interfaces 
such as those with atomic steps, valley-orbit mixing 
occur J^^^ ^^. Since both of the limiting cases exhibit 
strong dephasing due to intervalley phonons (with orbital 
excitations suppressed by a factor of 8 from valley excita- 
tions; see Eqs. ([ 



mixing cannot 
phasing. 



32) and (33)), we expect that valley-orbit 



3e used to suppress substantially this de- 



IV. DEPHASING DUE TO CHARGE NOISE 

We now consider the other expected major dephasing 
mechanism for our system: charge noise. ^^^^^^1^ This de- 
phasing arises because remote charge fluctuations induce 
random variations of the energy splitting between sin- 
glet and triplet levels by coupling to their non-equivalent 
charge distributions via the Coulomb interaction. These 
variations in the energy splitting lead to the accumula- 
tion of a random phase between the singlet and triplet 
states. In turn, this introduces a phase difference be- 
tween the logical qubit states. The dephasing mechanism 
and the estimated decay rates are essentially equivalent 
for GaAs and Si quantum dots, so we do not treat these 
two systems separately in this section. 

Because the hybrid qubit has a potentially strong 
charge characteristic, it can couple to remote charge 
traps. We assume the simplest no n- trivial charge fluc- 
tuation: a single, remote charge trap with two states 
(occupied and empty). To determine the dephasing rate, 
the first step is to compute the effect of the change in 
the state of the charge trap on the singlet-triplet energy 
splitting. We work to first order in perturbation theory, 
where we may calculate the small change in energy by us- 
ing the unperturbed (spatial) wavefunctions. Using the 
formalism of Sec. [llj it is straightforward to show that 
the first-order estimate of the variation in energy split- 
ting AV{r) is 



AV{r) = {T\V{r)\T) 
«(l|U(r)|l)- 



{S\V{r)\S) 
;0|U(r)|0), (34) 



where |0) is the ground state, |1) is the first excited state, 
and r is time. Here, we assume that the energy fluctua- 
tions are much smaller than the singlet-triplet splitting. 



and hence also the confinement energy. At any instant 
in time, our perturbing charge trap might be occupied or 
empty. If the charge trap at the position r is occupied, 
and hence perturbing the singlet-triplet energy splitting, 
we have 



AU: 



/.3. 



|^i(r')r-|^o(r')r 



47reo 



(35) 



where e is the elementary charge and eo is the (low- 
frequency) dielectric constant of our material. We as- 
sume that the trap is distant and calculate AU in a multi- 
pole expansion. ^^ For S- and P-like orbitals in a perfectly 
harmonic dot, the lowest-order, non- vanishing term is of 
quadrupole order: 



AFsp 



32 l2 (1 
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(36) 



where L is the lateral electron confinement length and 
(r, ^, (j)) is the location of the noise source in polar- 
spherical coordinates. Alternatively, for two valley states 

with identical envelope functions, we find that AU 

,2 7 2 
is exponentially suppressed by a factor of e~^ '^o to 

quadrupole order, where d is the z-envelope width and 
ko is the location of the valley minimum. To evaluate 
Eq. (36), we must estimate the typical distance r be- 
tween the charge trap and the qubit. To do this, we 
consider a slightly different system comprised of a dou- 
ble dot charge qubit in GaAs, for which the energy split- 
ting has been measured experimentally and found to be 
^yexp ^ I Q ^eV.^^ Although this system differs from a 
two-electron dot, the statistics of the charge fiuctuators 
should be similar. Evaluating Eq. (35) for the double-dot 



geometry, we find that the leading order term is a dipole 
contribution: 



AU; 
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(37) 



where po is the dipole moment exo associated with 
the dot separation xq. Averaging over 6 and (j) with 
xq ~ 300 nm (the distance between the double dots con- 
sidered in the experiment), we solve Eq. (37) to find 
r ~ 2.9 /im. Inserting this into Eq. (36), we obtain 
AU^p ^2x 10"^ /ieV for L = 40 nm in Si. Note that this 
value for AU^-p is likely an overestimate: if the dephas- 
ing is due to multiple charge traps (instead of the single 
trap we have assumed), AU^-p will be decreased. This is 
because matching to AU^^^ while increasing the number 
of traps increases the average r. Since AUs-p ~ 1/r^ and 
^Vexp ^ l/^^7 AU^-p decreases. Therefore, our estimate 
of the S — T dephasing rate for an ideal orbital first ex- 
cited state is an overestimate, and may decrease due to 
multiple charge traps. 

Now that we know the magnitude of the energy fluctu- 
ations, we can calculate the dephasing time T2. The off- 
diagonal elements of the density matrix decay as e~^^^^\ 
so the time T2 is defined by A(/)(T2) = 1.^^ Following 



Ref. [171 the time-dependent dephasing is given by 



A<j>iT) = ^ £ d^Sico) 



smujr/2 
uj/2 



(38) 



where cjq is a low- frequency cutoff that is the inverse mea- 
surement time. Up to this point, we have only considered 
the coherent evolution of the phase due to charge fluc- 
tuations. However, true decoherence occurs due to the 
statistical nature of the fluctuators. This effect is cap- 
tured in the spectral density S{uj) of the charge noise, 
through the definition 



1 r 



dre'^^ (AV(r)AV(O)) . (39) 



As noted in Ref. 17^ by examining the form of Eqs. (38) 



and (39) we can deduce that T2 for 1// noise should 
scale as 1/AV. Thus, we can calibrate our T2 to the 
experimental measurement via 
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where A V^^^ and T'2 ^^ are the experimental charge qubit 
measurements for the energy splitting fluctuation and 
the dephasing time. For the double-dot charge qubit ex- 
periment referenced above, T^^^ ?^ 1 ns, which leads to 
T2^ ^ 0.8 /is for our two electron dot with an orbital- 
like first excited state. Thus, the dephasing rate for or- 
bital excitations due to charge noise in Si is F^^^^^^ ^ 
1.3 MHz, which is on the same order as the phonon- 
induced dephasing in Si, and far slower than the dom- 
inant dephasing mechanism in GaAs. We plot this de- 
phasing rate as a function of L alongside the phonon 
dephasing mechanisms in Fig. [l] 

This long dephasing time is due to the fact that for 
a perfectly harmonic confinement potential, the dipole 
term of AVsp vanishes. However, in realistic systems, 
potential anharmonicity and interface roughness result 
in a non- vanishing dipole moment that can be more im- 
portant than the quadrupole term in Eq. (36). If our 



confinement potential is severely anharmonic, we expect 
that we would have a dipole contribution similar to Eq. 
(37), but with a moment of po ^ eL. Further, a dis- 



ordered interface might also introduce a dipole moment. 
As an example, simulations of the typical devices used in 
Refs.|46]and|47|findpo/e = 1.8 nm. 

Fig.[3]shows the dephasing rate due to charge noise as a 
function of dipole moment po, obtained by using Eqs. ( 37 ) 
and (pOl). The figure also shows the dominant dephasing 



rates from phonons calculated in Sec. Ill and listed in 
Table [Tj which are essentially independent of dipole mo- 
ment because there is substantial electron-phonon cou- 
pling even for perfectly harmonic confinement poten- 
tials. Within our approximations, we see that phonon- 
mediated dephasing is the most important mechanism in 
GaAs, but in Si charge noise can easily dominate. 



FIG. 3: (Color online) Singlet-triplet dephasing rate Tst 
due to charge noise and electron-phonon coupling as a func- 
tion of effective dipole moment po- In this plot, a constant 
quadrupole contribution of 1.3 MHz (Eq. (36)) is added to 
the dipole contribution, which is the estimated dephasing rate 
from charge noise for a dot with a purely harmonic confine- 
ment potential, for which po/2 is zero. Here, we have set the 
lateral electron confinement to be L = 40 nm, the vertical con- 
finement to be cZ = 3 nm, and have assumed an orbital excited 
state. The charge noise curve for Si is estimated using a di- 
electric constant eg* = 11.7e^ac. In a perfectly harmonic dot, 
po ~ 0, but anharmonicity and disorder can introduce a dipole 
moment of po ^ eL. In GaAs we expect phonon-mediated de- 
phasing to be most important, but in Si quantum dots charge 
noise can easily dominate. The circle marker indicates the 
estimated dephasing due to charge noise at vo/ e — 1 .8 nm, 
an estimated dipole moment for realistic devices .'^^'^ 



V. DISCUSSION 

In this paper, we addressed dephasing due to electron- 
phonon coupling and charge noise for two-electron states 
in a single quantum dot in both GaAs and Si. For the 
electron-phonon coupling, we found that in GaAs the 
main contribution to dephasing is due to polar coupling 
to optical phonons, and that the dephasing rate was of 
order gigahertz. In Si, phonon-mediated dephasing rates 
are much lower than in GaAs because there is no po- 
lar coupling to phonons, Intervalley processes are more 
important than intravalley processes, since phonons that 
couple valleys in silicon are extremely short-lived. The 
intervalley coupling to phonons leads a dephasing rate 
for silicon of order megahertz. 

We found that charge noise for an orbital first excited 



state in a perfectly harmonic quantum dot with no dis- 
order is strongly suppressed because the singlet-triplet 
energy splitting fluctuations produced by a remote per- 
turbing potential in this case are of quadrupole order, 
while for a double-dot charge qubit they are of dipole or- 
der. As has been noted previously,^^ pure valley states 
in Si are even more favorable, as they are largely immune 
to charge noise, in that both the dipole and quadrupole 
terms are suppressed by a factor of e~^ ^o. However, the 
introduction of either anharmonicity (for orbital excita- 
tions) or disorder (for valley excitations) leads to non- 
vanishing dipole moments up to the order of the lateral 
wavefunction extent. For either type of excited state in 
Si, this can become the dominant dephasing mechanism. 
Our estimated dephasing rate due to charge noise, based 
on calculations in typical dots, is of order 10 MHz. This 
rate is fast enough to dominate the dephasing in silicon, 
but likely not in GaAs. 

Our calculations suggest that two-electron, singlet-dot 
systems in Si have substantially better dephasing prop- 
erties than those in GaAs. This is because the polar cou- 
pling for optical phonons, which mediates fast dephasing 
in GaAs, is absent in Si. Within Si, to reduce the dephas- 



ing in this system, the critical figure to optimize is the 
effective dipole moment of the charge density difference 
between the first excited state and the ground state. As 
indicated in Fig. [3J we estimate that unless this dipole 
moment is reduced below po/e ^ 1 nm, charge noise is ex- 
pected to be the dominant dephasing mechanism. Below 
that threshold, the electron-phonon coupling (for valley 
excitations) and quadrupole-order charge noise (for or- 
bital excitations) are expected to be the dominant de- 
phasing mechanisms. 
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